Excitation spectrum and effective interactions of highly-elongated Fermi gas 
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Full 3D calculations of small two-component Fermi gases under highly-elongated confinement, 
in which unlike fermions interact through short-range potentials with variable atom-atom s-wave 
scattering length, are performed using the correlated Gaussian approach. In addition, microscopic 
ID calculations are performed for effective "atomic" and "molecular" ID model Hamiltonian. Com- 
parisons of the 3D and ID energies and excitation frequencies establish the validity regimes of the 
effective ID Hamiltonian. Our numerical results for three- and four-particle systems suggest that 
the effective ID atom-dimer and dimer-dimer interactions are to a good approximation determined 
by simple analytical expressions. Implications for the description of quasi-lD Fermi gases within 
strict ID frameworks are discussed. 
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I. INTRODUCTION 

Ultracold atomic and molecular gases arc considered 
nearly ideal model systems since their confining geome- 
try, size and interaction strength can be varied with un- 
precedented control 0]. A key goal of ongoing research 
activities is to experimentally determine the complete 
phase diagram of cold atom systems 01 . The success- 
ful demonstration of this task would provide a first step 
towards utilizing cold atom systems as quantum emula- 
tors. The determination of phase diagrams of quasi-lD 
systems has received considerable attention since these 
systems can, under certain circumstances, be described 
by ID model Hamiltonian whose properties have been 
studied extensively in the literature Q. For this class 
of systems, the challange is to establish which aspects of 
quasi-lD cold atom experiments can be described by ID 
model Hamiltonian. 

Naively, the ID scattering strength between two par- 
ticles in a wave guide geometry may be estimated by 
integrating out the tightly-confined transverse degrees of 
freedom. However, while the result is accurate in the 
weakly-interacting regime, Olshanii's seminal work [4| 
shows that the ID scattering strength gf^ depends in 
general non-trivially on the 3D s-wave atom-atom scat- 
tering length and the transverse angular frequency 
Up. The coupling constant gffy determined by Olshanii is 
now widely used in many-body studies of Bose and Fermi 
gases. The applicability of an effective atomic ID Hamil- 
tonian whose two-body interactions are parameterized in 
terms of has, e.g., been confirmed for a Bose gas un- 
der highly elongated harmonic confinement by comparing 
the results of 3D and ID Monte Carlo calculations [5[. 

Over the past few years, an effective atomic ID 
Hamiltonian has also been applied extensively to two- 
component Fermi gases under highly-elongated confine- 
ment 0, 0, in this case, however, the validity regime 
of the effective atomic ID Hamiltonian has not yet been 
assessed carefully. It is clear that an effective atomic 
ID Hamiltonian description breaks down when tightly- 
bound molecules form. In this case, the system may be 



described by an effective molecular ID Hamiltonian that 
treats each tightly-bound molecule as a composite boson. 
While the functional form of such an effective molecular 
Hamiltonian is generally agreed upon, the parametriza- 
tion of the effective atom-dimer and dimer-dimer interac- 
tions varies 0,0]. Furthermore, it is not clear whether or 
not the effective atomic and molecular ID Hamiltonian 
descriptions connect smoothly in the strongly-interacting 
regime. 

This work presents 3D and ID zero-temperature ab 
initio calculations for small two-component Fermi gases 
with up to N — 4 atoms under highly-elongated confine- 
ment and assesses the validity regimes of effective atomic 
and molecular ID Hamiltonian. Our main findings are: 
i) The 3D energies are reproduced well by an effective 
atomic ID Hamiltonian for small \a^/a p \ (a|J, < 0), 



where a p denotes the oscillator length in the tight con- 
finement direction [see Eq. flSJ)]. ii) For small positive 
a 3E)/ a p> tne 3D energies are reproduced well by an effec- 
tive molecular ID Hamiltonian that depends on the effec- 
tive ID atom-dimer and dimer-dimer scattering lengths 
afp and a^; analytical expressions for and are 
presented. Hi) For two of the energy curves considered 
(see below) , the descriptions based on the effective atomic 
and molecular ID Hamiltonian join fairly smoothly in the 
strongly-interacting regime, defined through \a^\ > a p ; 
not surprisingly, the dependence of the energies on the 
aspect ratio is largest in the strongly-interacting regime. 

Our assessment of the validity regimes of the effective 
atomic and molecular ID Hamiltonian for small systems 
is expected to provide guidelines for larger systems, and 
is thus of great importance for realizing condensed matter 
and materials analogs as well as for exploiting cold atom 
systems for quantum computation and quantum simu- 
lation. Quasi-lD few-fermion systems can be prepared 
by loading a gas of ultracold fermions into an optical 
lattice [9[. Measurements of the excitation spectrum as 
a function of the interaction strength would provide a 
stringent test of our microscopic predictions. 

Section [Hi introduces the 3D model Hamiltonian, dis- 
cusses the numerical techniques employed to solve the 
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corresponding Schrodinger equation and presents the re- 
sulting 3D energies. Section HU introduces the effective 
atomic and molecular ID Hamiltonian and presents de- 
tailed comparisons between the 3D and ID energies. Sec- 
tion IIVI discusses the excitation spectrum of strongly- 
interacting two-component Fermi gases under highly- 
elongated cylindrically-symmetric confinement. Finally, 
Sec. IVl concludes. 



II. FULL 3D TREATMENT: ENERGETICS 

This section introduces the 3D model Hamiltonian and 
the numerical techniques employed to solve the corre- 
sponding Schrodinger equation. 3D energies are pre- 
sented for N = 2 — 4 fermions under highly-elongated 
confinement. 

Our 3D model Hamiltonian H 3 d for the trapped two- 
component Fermi gas with N\ spin-up and N2 spin-down 
fermions, where N = N\ + N2, reads 
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Here, m and fj denote the atom mass and the position 
vector of the ith atom, fj = (xi,yi,Zi). The trapping 
potential Vtr(^) is given by 
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(2) 



where pi and A are defined through pi — \J x \ + yf and 



(3) 



and ui p and to z denote the transverse and axial angular 
frequencies. Unlike atoms interact through a spherically 
symmetric short-range Gaussian potential Vtb) 



Vtbinj) =-y exp ( 



(4) 



where = \ri — fj\. We take the range ro to be much 
smaller than the oscillator lengths a z and a p in the z- 
and p-directions, 



2 ,p 



(5) 



and adjust the depth Vq (Vq > 0) so that the free-space 
3D s-wave atom-atom scattering length a^y takes the de- 
sired value. A solid line in Fig. [1] shows function 
of the well depth Vq. To realize different negative Ogp, 
we start with a non-interacting (NI) system (Vq = 0) and 
increase the depth Vb till |a^| becomes infinitely large; 
at this point, the free-space two-particle system supports 
a single zero-energy s-wave bound state. To realize dif- 
ferent positive dgp, we increase Vq further. In general, 
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FIG. 1: (Color online) The solid and dashed lines show the 
free-space scattering lengths and a aa ) for the Gaussian 
potential Vtb as a function of the well depth Vb. The scattering 
lengths and depths are measured in the "natural units" of 
the free-space system, i.e., in units of ro [zo] and /i 2 /(2/i aa rQ) 
{h 2 /(2p^zl)\ for the 3D [ID] system, where ^ denotes the 
reduced mass of the atom-atom system. 



Vtb can lead not only to s-wave scattering but also to 
higher partial wave scattering. We have checked that 
the generalized p-wave scattering length and generalized 
scattering lengths corresponding to other higher partial 
waves are negligible over the range of well depths con- 
sidered in this paper. This implies that H^> effectively 
describes an s-wave interacting system. 

To solve the time-independent Schrodinger equation 
for i?3D, we separate off the center-of-mass motion and 
numerically solve the resulting Schrodinger equation in 
the relative coordinates. For N\ = N2 = 1, we expand 
the relative wave function in terms of two-dimensional B- 
splines and diagonalize the Hamiltonian matrix. For the 
three- and four-particle systems, we employ a correlated 
Gaussian (CG) approach [l(| EH that expands the rela- 
tive wave function ip in terms of Gaussian basis functions 
Ap) Az) 

Jk Jk ' 

Nb 



j> = ^ CkA 
fe=i 

/fc P) (Pl2, ' ' ' ,PN-l,N)f ( k Z \zi2,--- ,ZN-1 x n)] , (6) 



where 

/fc P) (pi2, • • • ,Pn-i,n) = exp 



Pij 



N 

— I x/2d {p) 

K] \ VZ,a ij,k, 



(7) 



(z) 

and fj- ' is defined analogously. The relative coordinates 
Pij and Zij are defined as pij = \J (xi — Xj) 2 + (jji — yj) 2 
and Zij = Zi - 

d) p \ and d<*\ are chosen semi-stochastically for each pair 



■Zj = !,••• ,N with i < j). The widths 



ij,k ij^k 

ij and fcth basis function, and the total number of basis 
functions is denoted by Nb- In Eq. ([6]), the Ck denote 
expansion coefficients and A denotes an anti-symmetrizer 
that ensures the proper symmetry of the two-component 
Fermi gas under exchange of identical fermions. For N = 



3 (Ni = 2 and N2 = 1), A can be conveniently written 
as 1 — P12, where Pyi permutes the two up-fcrmions. For 
N = 4 (Ni = N 2 = 2), A can be written as 1 - P i2 - 

For the interaction and confining potentials chosen, the 
Hamiltonian and overlap matrix elements (the basis func- 
tions fj, fjf^ do not form an orthogonal set) can be con- 
structed analytically. The diagonalization of the eigen- 
value equation is then performed using standard tech- 
niques. The resulting eigenenergies, whose accuracy can 
be systematically improved by increasing the number of 
basis functions and by optimizing the widths d$ k and 

(z) 

"'ij k °f ^ ne Gaussian functions, provide upper bounds to 
the exact eigenenergies. 

The 2D functions fjf ] are eigenfunctions of the z- 
component L z of the orbital angular momentum operator 
with eigenvalues hmi , mi — [101 ] , while the ID functions 

(z) 

ft have even parity P z = +1. For the N = 4 sys- 
tem, the energetically lowest-lying state has to/ = and 
P z = +1 for all 3D scattering lengths ag, and the basis 

functions fjf fi defined in and below Eq. ((6]) have the 
proper symmetry. The ground state of the NI N = 3 sys- 
tem, in contrast, has to./ = and odd parity (P z = —1), 

which cannot be described by the basis functions fff* f z z ^ . 
To describe states with odd parity, we add a spectator 
atom that does not interact with the ./V-fermion system 
of interest; the energy of the NI spectator atom follows 
from its to/ quantum number and from its parity, and is 
subtracted at the end of the calculation. Since the basis 
functions of the (7V+l)-system have even parity, the spec- 
tator atom and the 7V-fermion system either both have 
even parity or both have odd parity. In the following, we 
label our solutions by the parity P z ; if a spectator atom is 
added for computational purposes, we report the parity 
of the physical system of interest. Furthermore, since all 
energetically lowest-lying states of two-component Fermi 
gases under highly-elongated confinement have m/ = 0, 
we frequently omit the to; label. 

Figures (Ha) and (b) show the relative 3D energies for 
N = 2 — 4 and A = 10 as a function of the inverse 3D 
scattering length a p /a^y. Solid lines in Figs.^a) and (b) 
show the relative s-wave ground state energy E s (l,l) — 
Ecm,o of the two-body system. For later convenience, the 
energy £7,(1,1) as well as the energies E(N\,N2) (see 
below) include the center-of-mass ground state energy 
Ecm,o, E C mm = hu p + hiu z /2. In the NI limit, E„(l, 1) — 
£cm,o equals 10.5huj z = l.Qbhwp. In the absence of the 
confining potential in the z-direction, the relative two- 
body energy is always smaller than Hu> pi indicating the 
existence of a quasi-lD bound state for all 3D scattering 
lengths [J [3| ■ The confining potential in the in- 
direction pushes the energy up; in the NI limit, the up- 
shift is given by the zero-point energy Huj z /2. 

For the N = 3 system, the relative energies E(2, 1) — 
Ecm,o of the energetically lowest-lying states with P z = 
+1 and —1 are shown by dotted and dashed lines, re- 
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FIG. 2: (Color online) Relative 3D energies calculated us- 
ing Hzu as a function of dp/a^ for (a) negative and 
(b) positive ag, for A = 10 [r = 0.03a z for ag, < and 
N = 4, and 0.01a z otherwise]. Solid lines show the s-wave en- 
ergy E s (l, 1) — -Ecm.o of the two-particle system, dashed lines 
show the lowest three-particle energy E(2, 1) — Ecm,q with 
P z = — 1, dotted lines show the lowest three-particle energy 
E(2, 1) — Ecm,o with P z — +1 and dash-dotted lines show 
one half of the lowest four-particle energy E(2,2) — Ecm,o 
with P z = +1. The horizontal solid lines on the left side 
of panel (a) indicate the relative energies of the NI sys- 
tems: E s (l, 1) - Ecm,o = W.5hiu z , E(2, 1) - E CM ,o = 22}uj z 
(P z = -1), £(2,1) - Ecm.o = 23hujz (Pz = +1), and 
[E(2, 2) - £cm,o]/2 = 16.75hwz. 



spectively. The P z — — 1 state has lower energy for small 



'3D I 



a 3D < [see Fig. Ufa)], while the P z = +1 state 



has lower energy for small positive 



'3D 



[the 



crossover 



of the two states is not visible on the scale shown in 



Fig. [2b); it occurs at a p /a 



aa 
3D 



2 (see also Fig. 0J] 



The relative three-particle energies are just slightly larger 
than the relative two-body s-wave energies in the limit of 
small positive , indicating that the three-particle sys- 
tem can be thought of as consisting of an s-wavc dimcr 
and an unpaired atom. The relative ground state energy 
E(2, 2) — Ecm.o of the four-particle system has P z = +1 
for all 3D scattering lengths ; to ease comparisons be- 
tween the energies of the two- and four-particle systems, 
dash-dotted lines in Figs. E^a) and (b) show one half of 
the relative four-particle energy. For small positive a^, 
the four-particle energy approaches approximately twice 
the energy of the two-particle system, indicating that 
the four-particle system can be thought of as consist- 
ing of two s-wave molecules. No tightly-bound trimers 
or tetramers are formed in the — > + limit, in agree- 
ment with results for zero-range interactions [H, [lrlll7j . 

The 3D energies can be combined to define the univer- 
sal energy curve Ajv^jVa for a two-component Fermi gas 
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FIG. 3: (Color online) Ajvi.jVa as a function of ap/a^: 
Panel (a) covers a large scattering length range while panel (b) 
shows an enlargement of the strongly-interacting regime. 
Thick dotted, dashed and solid lines show A2,2 (Pz = +1), 
A2,i (Pz = +1) and Aa,i (P z = —1) calculated using H30 for 
A = 10 [r = 0.03a z for a|f) < and N = 4, and 0.01a z other- 
wise]. Thin dash-dash-dotted and dash-dot-dotted lines show 
the corresponding ID energy curves calculated using Hf u and 

rrm 
#1D- 



under external cylindrically symmetric confinement, 



A 



Nt,N 2 



E(N U N 2 ) - NdEJX 1) - Nfhjwp + \w z ) 
E m (N 1 ,N 2 )-Nh(w p + ^u; z ) 



where N d = min{N 1 ,N 2 } and N{ = \Ni~N 2 \. In Eq. ©, 
Em (Ni , iVa) denotes the energy of the NI system, and the 
energies Em(Ni,N 2 ), E(Ni,N 2 ) and E s (l, 1) include the 
center-of-mass ground state energy Ecm,0- To remove de- 
pendencies of the total energy E(Ni,N 2 ) of the trapped 
system on 14b, the s-wave ground state energy E s (l,l) of 
the trapped two-particle system is subtracted on the right 
hand side of Eq. ([5]). If E(Ni,N 2 ) corresponds to the en- 
ergetically lowest- lying state of the NI system, the univer- 
sal energy curve Ajv x ,n 2 equals one. Conversely, if A.ff u N 3 
equals zero in the — > + limit, then this indicates 
that the system is effectively NI and that induced inter- 
actions are absent. The definition of the universal energy 
curve presented in Eq. (fSJ) for cylindrically-symmetric 
two-component Fermi systems constitutes a straightfor- 
ward generalization of that previously introduced for 
spherically-symmetric two-component systems [H, HH • 

Figure [3] shows the 3D energy curves An 1} n 2 as a func- 
tion of cip/a^ for A = 10; the energies used to calculate 
the AjVi,jv a are the same as those shown in Fig. [51 A thick 
dotted line shows A 2y2 calculated using the four-body en- 



ergies E(2, 2) that correspond to states with mi = and 
P z = +1. The energy curve h. 2 ^ 2 decreases monotoni- 
cally from 1 to approximately as ap/a^ increases from 
small negative to large positive values. Thick dashed and 
solid lines in Fig. [3] show the energy curves K 2 ^i for the 
energetically lowest-lying N = 3 states with P z = +1 
and —1, respectively. It can be seen that these N = 3 
energy curves cross at a p /a|fj « 1.9. For — > (NI 
limit), the ground state has mi = and P z = — 1: One 
spin-up and one spin-down atom occupy the ground state 
harmonic oscillator orbital while the second up-atom oc- 
cupies the first excited state harmonic oscillator orbital. 
For a|f) — > + , in contrast, the ground state for N = 3 
has mi — and P z = +1: The system consists of a 
tightly-bound dimer and an atom, which both occupy 
the lowest trap state. In this limit, the energy of the 
state with P z = — 1 is about huj z larger than the energy 
of the P z = +1 state [note that the energy difference is 
too small to be visible on the scale shown in Fig. G^b)]. 
This suggests that the tightly-bound molecule and the 
unpaired atom interact through effective ID potentials 
that lead to even and odd parity scattering for P z = +1 
and — 1, respectively (see also the next section). 



III. ID TREATMENT: ENERGETICS AND 
EFFECTIVE INTERACTIONS 

This section considers effective atomic and molecular 
ID Hamiltonian, which assume that the motion in the 
p-direction is frozen, and compares the resulting ID en- 
ergies with the 3D energies discussed in the previous sec- 
tion. The applicability of the ID model Hamiltonian and 
their parametrizations are discussed in detail. 

If the system behaves like an atomic gas, the effective 



(8) atomic ID Hamiltonian is given by 



H 



ID 



where 



N r 

E 

i=i l 



-2m-dz}+ VtM) 



E 

1=1 j: 



N 

E 

=iVi + l 



Vtr (zi) = -muj 2 z zf. 



(10) 



In Eq. @, the spin- up and spin-down fermions interact 
through the two-body potential Vtbizij) and, as in the 
3D Hamiltonian [see Eq. (pQ)], like atoms do not 

interact. The two-body potential 14b is chosen such that 
its ID even parity atom-atom scattering length is 
given by [|| 



1 



a lD = 



2^ a 



- 1.4603 



l 3D 



(11) 



where /i aa denotes the reduced mass of the atom-atom 
system. The ID even parity scattering length a a f) [solid 
lines in Figs. Ufa) and (b)] is large for large lap/agpl 
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FIG. 4: (Color online) ID even parity atom-atom scattering 
length affj/ttp (solid line) and corresponding ID scattering 
strength gffy/(frw p a p ) (dashed line) as a function of the in- 
verse 3D atom-atom scattering length a p /a^fy: Panel (a) cov- 
ers a large 3D scattering length range while panel (b) shows 
an enlargement of the strongly-interacting regime. The ID 
coupling constant gffy changes from about —hui p a p to about 
— 60hu p a p as a p /a^fy increases from —1 to 1 (i.e., in the 
strongly- interacting regime). 

( a 3D < 0)i decreases linearly with increasing a p /a^fy, and 
crosses zero at a p /a^fy = 1.0326. Since the ID coupling 
constant g^, 

h 2 

rf* = , (12) 

diverges when vanishes, the quasi-lD system is in- 
finitely strongly- interacting for a finite a^fy. Further- 
more, a large positive a^fy indicates the presence of a 
weakly-bound even parity two-body bound state. In the 
literature, the effective ID atom-atom potential Vtb is 
frequently modeled by a ID zero-range (5-function poten- 
tial. For numerical convenience, we use instead a ID 
Gaussian potential [Eq. Q with and ro replaced by 
Zij and zq] with a small width zq (zq — 0.005a z ) and a 
depth adjusted so as to obtain the desired a^fy. We have 
checked that the resulting ID energies depend only very 
weakly on zq and that the odd parity atom-atom scat- 
tering length is negligibly small over the range of well 
depths considered. The ID even parity atom-atom scat- 
tering length a^p for the ID Gaussian potential is shown 
in Fig. [1] by a dashed line as a function of the well depth 

v . 

For the 3D energy curves considered in Fig. [31 the ef- 
fective atomic ID Hamiltonian is expected to provide an 
accurate description if the size of the ID dimer is much 



larger than the oscillator length a p (see, e.g., Ref. [1]). 
Approximating the size of the dimer by |a^| and using 
l a iE)l ~ l a p/ a 3E>l [i- e -; using the first part on the right 
hand side of Eq. (fTTj) ]. the validity condition a p ^> \a^\ 
is obtained. Relaxing the disparity of length scales, we 
have a p > logfjl with ag, < 0. 

To obtain the ID energies of the effective atomic ID 
Hamiltonian .fff D , we first separate off the center-of-mass 
motion and then solve the resulting Schrodinger equation 
in the relative coordinates using the B-spline approach 
for the two-particle system and the CG approach for the 
three- and four-particle systems. Our CG implementa- 
tion for the ID system parallels that discussed above for 
the 3D system. The main difference is that the basis 
functions are now given by fjf ^ instead of by fjf^fjf^- 
Systems with odd parity are, similarly to the 3D case, 
treated by adding a NI spectator atom. For large positive 
a^p, we find that the energetically lowest-lying ID states 
accurately model the corresponding 3D states. For small 
positive a^, however, the effective atomic ID Hamil- 
tonian supports a sequence of tightly-bound three- and 
four-particle states, which have no analog in the 3D sys- 
tem (as discussed above, tightly-bound three- and four- 
particle states are not supported by ^d); these ID states 
are excluded from our analysis. The ID energy states 
of interest to us are those that smoothly evolve from a 
NI gas-like state in the a^fy — > oo limit to states that 
describe a weakly-bound molecule and an atom or two 
weakly-bound molecules for N = 3 and N = 4, respec- 
tively, in the a^fy — > + limit. Our last ID energies are 
reported for a p /a^fy w 1 for N — 3 with P z — +1, and 
for ap/ag) = for N = 3 with P z = -1 and N = 4 with 
P z = +1. We note that tightly-bound ./V-body states 
also exist for negative . Their existence can be traced 
back to the finite range of the Gaussian two-body interac- 
tion potential; an effective atomic ID Hamiltonian with 
zero-range (5-function potentials and negative a^fy does 
not support tightly-bound TV-body states. 

The ID energies determine the ID energy curves 
Ajv^/Vj, which are given by Eq. (|SJ) with u p = 0, E s (l, 1) 
denoting the ID even parity two-particle energy, and 
E(Ni,N 2 ) and Em{Ni,N2) interpreted as ID energies. 
Thin dash-dash-dotted lines in Fig. [3] show the ID energy 
curves obtained using Hf D for N = 3 and 4. The agree- 
ment between the ID energy curves and the correspond- 
ing 3D quantities (thick lines) in the weakly-attractive 
regime is excellent. For P z = +1, the agreement be- 
tween the ID and 3D energy curves extends into the 
strongly-interacting, large positive a^fy regime. The ID 
energy curve for N = 3 with P z = — 1, in contrast, starts 
deviating from the corresponding 3D energy curve for 
somewhat less strong interactions (\ap/a^,\ < 2 with 
a 3D < 0). We have checked that these deviations are 
not due to the finite range of I4t>- 

In addition to an effective atomic ID Hamiltonian H^ B , 
we consider an effective molecular ID Hamiltonian H^. 
We show in the following that the 3D energy curves with 
P z = +1 can be reproduced well for a^fy > by treating 
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the N = 3 and 4 systems as effective two-particie systems 
that consist of an atom and a tightly-bound molecule and 
of two tightly-bound molecules, respectively. To this end, 
we model the atom-dimer and dimer-dimer interactions 
through a (5-function potential. The effective two-particle 
ID Hamiltonian for the relative coordinate z then reads 

H ™ = + ^ V + 5 ' D<5(Z) ' (13) 

where j — ad and dd for the atom-dimer and dimer- 
dimer system, respectively, and where the ID coupling 
strengths and are related to the ID scattering 
lengths a^p and [Eq. (jT2J) with aa replaced by ad 
and dd, respectively]. We approximate the effective ID 
atom-dimer and dimer-dimer scattering lengths a^, and 
a iD ^ the right hand side of Eq. (jlip with superscripts 
aa replaced by ad and dd, respectively. The 3D atom- 
dimer and dimer-dimer scattering lengths and a^n, 
in turn, are approximated by their free-space values [12], 

Eil, 

o§£ = 1.18a£ (14) 

and mmm 

a& = 0.608ag>. (15) 

Physically, this implies that molecules are formed in 3D 
and that their effective 3D interactions with atoms and 
other molecules are renormalized by the quasi- ID con- 
finement. 

The validity regime of the effective molecular ID 
Hamiltonian Hf^ is expected to be determined by three 
conditions: i) Since the 3D free-space atom-dimer and 
dimer-dimer scattering lengths are derived assuming that 
a 3D r 0i the above parametrization is expected to 
break down when approaches r$. ii) For three- and 
four-particle systems under spherically symmetric con- 
finement, it has been shown [lj] that the full 3D energies 
on the BEC side (positive a^) are well described by 
effective 3D atom-molecule and molecule-molecule mod- 
els if eigp is much smaller than the harmonic oscillator 
length. Correspondingly, since our parametrization of 
the effective interactions given in Eqs. (fl"3")) - (fT5"|) for the 
highly-elongated system assumes that the molecules are 
formed in 3D, the validity regime of is expected 
to be given by <C a p . in) The effective ID model 
treats the dimer as a point particle. This treatment is 
justified if the atom-dimer and dimer-dimer distances are 
much larger than the size of the dimer, i.e., if a z ^> a^, 
(see, e.g., Ref. Q). Combining the three criteria, we find 
that -ffj^) is expected to provide an accurate description 
if a p 3> a^f) ^ ro or, employing less stringent criteria, 
if ftp ^ fl3D ^ ^o- Combining this with the expected va- 
lidity regime of (see above) , the strongly-interacting 
regime is defined through \a^\ > a p . If the exact effec- 
tive ID atom-dimer and dimer-dimer scattering lengths 
were known, condition ii) would not apply and the ex- 
pected validity regime of the molecular ID Hamiltonian 
would be larger (a z > a|p > r ). 



The Hamiltonian -fff^, given in Eq. (| 13[) parametrizes 
the effective interactions through a (5-function potential 
and thus assumes that the effective ID atom-dimer and 
dimer-dimer interactions lead to even parity scattering. 
Consequently, the Hamiltonian H"^ does not describe 
the P z = — 1 energy curve for N = 3. An effective molec- 
ular ID model for the N — 3 system with P z = — 1 
would include an effective ID interaction that leads to 
odd parit y sc attering such as a so-called zero-range 5'- 
potential [20]. Although interesting, an effective molec- 
ular ID description of the energy curve with P z = — 1 is 
not pursued in this work. 

Calculating the eigenenergies of H ^ from the known 
quantization condition [l9| , we find that the energy of the 
energetically lowest-lying state with gas-like character 
agrees well with the 3D quantities E(2, 1)—E S (1, l) — fhu p 
and £7(2,2) — 2E S (1, 1) for the atom-dimer and dimer- 
dimer systems with P z = +1 and > 0. Dash-dot- 
dotted lines in Fig.[3]show the energy curves for P z = +1 
calculated using the effective ID molecule model. These 
ID energy curves agree well with the corresponding 3D 
energy curves in the weakly-interacting regime. Devia- 
tions are visible for N = 3 for Op/a^ < 1 and for N = 4 
for a p /a^y < 0. The agreement of the ID and 3D en- 
ergy curves over a wide range of interaction strengths a 
posteriori justifies our parameterization of the effective 
ID atom-dimer and dimer-dimer scattering lengths (see 
also Sec. IIV[) . which differs from that employed in ear- 
lier work 0, ID, [TtJ • Notably, the ID energy curves for 
the effective molecular ID Hamiltonian connect nearly 
smoothly with those for the effective atomic ID Hamil- 
tonian in the strongly-interacting regime. 



IV. EXCITATION SPECTRUM 

This section discusses the behavior of the excitation 
frequency luq for systems with P z = +1. Within our 
3D framework, the excitation energy fiioo is defined as 
the difference between the first excited and the lowest 
P z = +1 states. The corresponding ID excitation energy 
is defined as the difference between the energies of the 
corresponding ID states. 

Circles in Fig. [5] show the 3D excitation frequency ojo 
for N = 3 with P z = +1 and A = 10. Panels (a) and 
(d) show ujq as a function of Op/a^,. The excitation 
frequency w equals 2uu z in the NI limit (a|£, — > 0~), 
reaches its maximum for infinitely large and its min- 
imum for Op/a^fy ~ 5, and increases monotonically to- 
wards 2oj z as l/ctgQ increases further. To illustrate the 
dependence of u>o on the aspect ratio, squares and dia- 
monds in Figs. [5jb) and (c) show the 3D excitation fre- 
quency luq for two larger aspect ratios, i.e., for A = 15 
and 20. Small dependencies of luq on A are visible in the 
strongly-interacting regime. 

To ease comparisons between the 3D excitation fre- 
quencies and those based on the ID Hamiltonian, 
Figs. E^b) and (c) show enlargements of the strongly- 
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5: (Color online) Excitation frequency ujq/uJz for N = 3 
P z = +1 and ro = 0.01a z . (a) u)o/u z as a function of 
10. Circles show the 3D excitation frequency, 



for A 



-"•3D 

while solid and dotted lines show the corresponding ID quan- 
tities calculated using iffp and H^, respectively. Panels (b), 
(c) and (d) show the strongly-interacting regime in more de- 
tail, (b) too/oJz as a function of — affj/fflz. Circles, squares and 
diamonds show the 3D excitation frequency for A = 10, 15 and 
20. A solid line shows the corresponding ID quantities calcu- 
lated using Hi D . (c) ujq/ijj z as a function of —a^,/a z . Circles, 
squares and diamonds show the 3D excitation frequency for 
A = 10, 15 and 20. A dotted line shows the corresponding ID 
quantities calculated using H^,. (d) Enlargement of panel 
(a). 



These scales are chosen since and determine the 



1 1D 



properties of and H^, respectively. Solid lines in 
Fig. [5] show the excitation frequencies calculated using 
the effective atomic ID Hamiltonian. These ID excita- 
tion frequencies reproduce the 3D excitation frequencies 
well in the weakly-attractive regime (a^ large). In the 
strongly-interacting regime, the agreement improves with 
increasing A [see Fig. [Sfb)]. Dotted lines in Fig. [5] show 
the ID excitation frequencies calculated using H^. For 
large |a^| (a^ < 0), the 3D excitation frequencies are 
independent of A and reproduced well by the ID molcc- 



*!dI 



the agreement improves 



ular model. For smaller 
with increasing A [see Fig. [5jc)] 

The effective molecular ID Hamiltonian predicts 
that a subset of the P z = +1 energy spectrum of the 
three- and four-particle systems coincides with that of 
a two-particle Tonks- Girardeau (TG) gas for a^> = 
[ij . For this atom-dimer scattering length, the effective 
molecular ID Hamiltonian predicts A2.1 = 1, A2,2 = 1/2 
and loq — 2co z . Assuming that the behavior of the effec- 
tive dimer system is indeed governed by (i.e., assum- 
ing that effective range and other corrections are negligi- 
ble), the condition a^, = signals an atom-dimer reso- 
nance. Our 3D calculations for the three-particle system 
with P z = +1 show that the ground state energy corre- 
sponds to that of a TG gas for a p /a^ w 1.54— 1.58 and 
that uj equals 2lo z for a p /a\v w 1.1—1.4 for A = 10 — 20, 
in fairly good agreement with the prediction based on 
the ID model, a p /a%v — 1.4069. The good agreement 
between our 3D results and those based on the effec- 
tive molecular ID Hamiltonian, which is based on a sim- 
ple empirical parametrization of the effective ID atom- 
dimer and dimer-dimer scattering lengths, suggests that 
the molecular ID model employed in this work provides 
a viable and fairly accurate description of the system. 

The atom-dimer s-wave resonance of quasi-lD systems 
found here, ap/a^ « 1.5, is somewhat smaller than that 
found by Mora et al. [ItJ by solving a set of inte gral equa- 
tions for zero-range interactions, dp/a^ ~ 1.85 [2l|. The 
3D Hamiltonian employed by Mora et al. accounts for 
the same physics as our 3D Hamiltonian and the deter- 
mination of the effective ID atom-dimer scattering length 
should, at least in principle, be exact [17J. It is not clear 
at present why our empirical molecular ID Hamiltonian 
provides a seemingly better description than Mora et aZ.'s 
results for » 0. 

We also analyzed the ground state energy and excita- 
tion spectrum for N = 4. The four-particle 3D energies 
are harder to converge than the three-particle 3D en- 
ergies, and comparisons between the full 3D excitation 
frequencies and the corresponding ID quantities are ac- 
companied by non-neglegible uncertainties. We find that 
our 3D results are consistent with the dimer-dimer s- 
wave resonance value predicted by [Eq. (fT3"|) with 

j = dd, and a^> given by Eq. (jlip with aa replaced by 
dd], ap/a^ = 0.89 [p. 
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V. CONCLUSIONS 

In summary, we have presented highly-accurate, micro- 
scopic 3D calculations for small highly-elongated Fermi 
gases with N = 2 — 4 atoms and reported the energies as a 
function of the interaction strength, covering the weakly- 
attractive and weakly-repulsive regimes as well as the 
strongly-interacting regime. In addition, the dependence 
of the energies on the aspect ratio was investigated for 
selected cases. While the role of the aspect ratio is negli- 
gible in the weakly-interacting regimes, its role becomes 
more important in the strongly-interacting regime, possi- 
bly indicating that vitual excitations of transverse modes 
become relevant. The full 3D energy curves with P z = +1 
are reproduced to a good approximation by effective 
atomic and molecular ID models whose effective inter- 
actions are given by simple analytical expressions that 
depend on the atom-atom s-wave scattering length a^, 
the aspect ratio A and the atom mass m. We find that the 
energies obtained from these effective atomic and molec- 
ular ID Hamiltonian join fairly smoothly in the strongly- 
interacting regime. Assuming that the effective ID atom- 



dimer and dimer-dimer scattering lengths govern the be- 
havior of the highly-elongated system, we deduced the 
positions of confinement-induced atom-dimer and dimer- 
dimer resonances from our energies. Whether the effec- 
tive ID models also connect fairly smoothly for larger 
systems is a pressing questions, in particular since the 
determination of the phase diagram of highly-elongated 
systems often times relies on strictly ID treatments. 

In the future, it will be interesting to extend the 
studies presented here to larger population-balanced 
and population-imbalanced two-component Fermi gases. 
While some microscopic calculations exist for strictly ID 
systems, microscopic 3D treatments that accurately ac- 
count for the dynamics along the tight and loose confin- 
ing directions are challenging. Furthermore, it will be 
interesting to compare the ID results obtained here for 
small systems with those obtained within the local den- 
sity approximation and to extend analogous comparisons 
to larger systems. 

Support by the NSF through grant PHY-0555316 is 
gratefully acknowledged. 
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